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Abstract 

The recently obtained solutions of the Dirac equation in the confining 
SU(3)-Yang-Mills field in Minkowski spacetime are applied to describe the 
energy spectrum of charmonium. The nonrelativistic limit is considered 
for the relativistic effects to be estimated in a self-consistent way and it 
is shown that the given effects could be extremely important for both the 
energy spectrum and the confinement mechanism. 
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1 Introduction 



Theory of quarkonium ranks high within hadron physics as the one of central 
sources of information about the quark interaction. Referring for more details 
to the recent up-to-date reviewilj, it should be noted here that at present some 
generally accepted relativistic model of quarkonium is absent. The description 
of quarkonium is actually implemented by nonrelativistic manner on the basis 
of the Schrodinger equation (concerning the general ideology here see, e. g., ref. 
U) and then one tries to include relativistic corrections in one or another way. 
Such an inclusion is not single-valued and varies in dependence of the point of 
view for different authors (see, e. g., ref.pi and references therein). It would 
be more consistent, to our mind, building a primordially relativistic model so 
that one can then pass on to the nonrelativistic one by the standard limiting 
transition and, thus, to estimate the relativistic effects in a self-consistent way. 
As follows from the main principles of quantum chromodynamics (QCD), the 
suitable relativistic models for a description of the relativistic bound states of 
quarkonium could consist in considering the solutions of the Dirac equation 
in a SU(3)-Yang-Mills field representing the gluonic field. Indeed, the Dirac 
equation in a SU (3)- Yang-Mills field is the direct consequence of the QCD La- 
grangian in the same way as the Dirac equation for the hydrogen atom is the 
direct consequence of the quantum electrodynamics (QED) Lagrangian. Fol- 
lowing the latter analogy, the mentioned SU(3)-Yang-Mills field should be the 
so-called confining solution of the corresponding Yang-Mills equations, i. e., it 
should model the quark confinement. Such solutions are usually supposed to 
contain at least one component of the mentioned SU(3)-field linear in r, the 
distance between quarks. Recently, in Ref. 4 a number of such solutions have 
been obtained and the corresponding spectrum of the Dirac equation describing 
the relativistic bound states in those confining SU(3)-Yang-Mills fields has been 
analysed. In this paper we should like to apply the results of ref. U to the 
description of the charmonium spectrum. Here we solve the inverse problem, i. 
e., we define the confining gluonic field components in the covariant description 
(SU(3)-connection) for charmonium (which corresponds to a potential of qq- 
interaction in the nonrelativistic description) employing the experimental data 
on the mentioned spectrumj^j. As a consequence, we shall not use any nonrela- 
tivistic potentials modelling confinement, for example, of the harmonic oscillator 
or funnel types, in particular, because the latter do not satisfy the Yang-Mills 
equations, while the SU(3)-gluonic field used by us does. Accordingly, in our 
case the approach is relativistic from the very outset and our considerations 
are essentially nonperturbative, since we shall not use any expansions in the 
coupling constant g or in any other parameters. 

Further we shall deal with the metric of the fiat Minkowski spacetime M that 
we write down (using the ordinary set of local spherical coordinates r, 'd, (p for 
the spatial part) in the form 

ds^ = g^.^dx^' ® dx" = dt^ - dr'^ - r^{di9^ + sin^ Mip^) . (1) 

Besides, we have \6\ — |det((7pi,)| = (r^sini?)^ and < r- < oo, < < tt, 
< < 27r. 

Throughout the paper we employ the system of units with h = c = 1, unless 
explicitly stated otherwise. Finally, we shall denote L2{F) the set of the modulo 
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square integrable complex functions on any manifold F furnished with an inte- 
gration measure while L2{F) will be the n-fold direct product of L2{F) endowed 
with the obvious scalar product. 



2 Preliminaries 

To formulate the results of ref.^ that we need here, let us notice that the 
relativistic wave function of the quarkonium can be chosen in the form tl) — 
{ipi, V'2)PSJ3) with the four-dimensional spinors ipj representing the j-th colour 
component of the quarkonium. The corresponding Dirac equation for ip may 
look as follows: 

= f^oi', (2) 

where jio is a mass parameter and one can consider it to be the reduced rela- 
tivistic mass which is equal, e. g. for quarkonia, to one half the current mass 
of quarks forming a quarkonium, while the coordinate r stands for the distance 
between quarks. 

From general considerations the explicit form of the operator T) in local coordi- 
nates x^^ on Minkowski manifold can be written as follows: 

V = ® h)E^{d^, ^h- ^cj^a67"7' ^h- igA^,), a < b, (3) 

where A — Af^dx^, A^ — A'j^Tc is a SU(3)-connection in the (trivial) three- 
dimensional bundle ^ over the Minkowski spacetime, I3 is the unit matrix 
3x3, the matrices Tc form a basis of the Lie algebra of SU(3) in the 3- 
dimensional space (we consider Ta to be Hermitian which is acceptable in 
physics), c = 1, 8, (g) here means tensorial product of matrices, g is a gauge 
coupling constant. At last the coefficients Ej^ and w^ab depend on the choice of 
metric and their form for metric (1) can be found in ref. 0] as well as the explicit 
presentation of matrices 7", a = 0,...,3. As for the connection A^ in bundle 
^, then, the suitable one should be the confining solution of the Yang-Mills 
equations 

d* F = *F A A - AA*F (4) 

with the exterior differential d — dtdt+drdr+d^)d'd+dipdip in coordinates t, r, ??, ip 
while the curvature matrix (field strentgh) for ^-bundle is = dA + A A A and 
* means the Hodge star operator conforming to metric (1). 
In ref.jlj the black hole physics techniques from refs.[n| were used to find a set 
of the confining solutions of eq. (4). For the aims of the given paper we need 
one of these solutions of ref.pP. Let us adduce it here putting Tc — Ac, where 
Ac are the Gell-Mann matrices (whose explicit form can be found in refs.'B'). 
Then the solution in question is the following one: 

A^+±AI = -^+A,,-A^+±Af = ^^-(A,+A,),--|a« = 

K+^K = bir+Bi,^Al+^Al -(6i+62)r-(Bi+S2), = ^2'^+^2 

(5) 

with all other = 0, where real constants aj, Aj, bj, Bj parametrize the solu- 
tion, and we wrote down the solution in the combinations that are just needed 
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to insert into (2). From the adduced form it is clear that the sohition is a con- 
figuration describing the electric Coulomb- like colour field (components At) and 
the magnetic colour field linear in r (components A^). Also, it is easy to check 
that the given solution satisfy the Lorentz gauge condition that can be written 
in the form div{A) = 0, where the divergence of the Lie algebra valued 1-form 
A = A^jTcdx^ is defined by the relation 



div(A) = -±^^d,{^\g^''A^) . (6) 

As was shown in ref.j3|, after inserting the above confining solution into Eq. 
(2), the latter admits the solutions of the form 

F,i(r)$,(^?,^) 



i/, = e'^^W-' ^'/V i'Va ^ ,J = 1,2,3 (7) 

with the 2D eigenspinor $j ~ ^ '^^ ^^"^ Euclidean Dirac operator on the 

unit sphere The explicit form of $j is not needed here and can be found in 
refs. |2] . For the purpose of the present paper it is sufficient to know that spinors 
<i>j can be subject to the normalization condition 

7T 2-JT 

I$ji|' + |$j2|')sini9di?d(p = l, (8) 



i. e., they form an orthonormal basis in L2(§^). 

The energy spectrum e of a quarkonium is given (in a more symmetrical form 
than in ref. 14j) by the relation e = uji + lu2 + lut, with 



-Aig^aibi + {ni + ai) v'(nf + 2mai + A?)^g + g'^hl[nl + 2niai) 
- ^i(ni,;i,Ai) nl + 2n,a,+Kl ' 

(9) 

UJ2 = t^2(«2,^2, A2) = 



-A2.g^(ai + a2)(bi + 62) - (^^2 + a2)\/{nl + 2n2a2 + ^^l)^Jl + g^jbi + b2Y(nl + 2n2a2) 

n1 + 2n2Q;2 + A^ 

(10) 

-A3g^a252 + (ng + a3)\/{nl + 2nzaz + A^)^g + g'^hljnl + 2^30:3) 



W3 = a;3(n3,;3, A3) 



n| -I- 2n3a3 -I- A§ 

(11) 

where Ai = Ai - gBi , A2 = A2 + ^(Bi + B2) ,A3 = A3 - 5S2 , Uj = 0,1,2,..., 
while Aj = i(^i + 1) arc the eigenvalues of euclidean Dirac operator on unit 
sphere with Ij = 0, 1, 2, ... Besides, 



ai = -y/A^ - 52^2 ^a2 = \lKl- g^{ai + 03)2 , 03 = y^A^ - ^2^2 _ (^2) 

Further, the radial part of (7), for instance, for the ■i/'i-component, is given at 
rti = by 



Fn - CiAr^^e-'^^'- [l-^\^F^2 = iC.Br"^ e-^^'' ( 1 + ^ ) , (13) 
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while at ni > by 



AB 



Fi2 = iCiSr"^e 



1 



r2ai/ \ j-2ai + l/ \ 

in) - (n) 



(14) 



with the Laguerre polynomials i^^(ri), ri = 2/3ir, /3i = \J [i^ — (t^i — g^i)^ + 
A = g6i+/3i, B = /xo+wi-gAi, = [ai/3i-gai(a)i-,gAi)+gai6i]B+g2Q(^^^^^ 

Z\ = [(Ai — 5-Bi)A + ga\\iQf\B + ^^ai^iA. Finally, C\ is determined from the 
normalization condition 



Jo 



(|Fn|' + |Fi2|')dr = -. 



(15) 



Analogous relations will hold true for ip2,3, respectively, by replacing ai, Ai,bi, Bi,ai 
a2,A2,b2, -82- 03 for ti'3 and ai, Ai,bi, B i, ai -(01 + 02), - {Ai + A2) , - {bi + 
b2),-{Bi+B2).02 ( '2 so that /32 = ^Mo - [^2 + g{Ai + A2)]^ + g'^{bi + 62)^ 



/Js = y'/ig — (cj3 — gA2Y + .9^fo|- Consequently, wc shall gain that -0^ G i|( 
at any t e M and, as a result, the solutions of (7) may describe relativistic bound 
states of a quarkonium with the energy spectrum (9)-(ll). 
Before apUying the above relations to a description of charmonium spectrum 
let us adduce the nonrelativistic limits (i e., at c ^ 00) for the energies of 
(9)-(ll). The common case is not needed to us in the present paper, so we shall 
restrict ourselves to the case of = 0, 1 and Ij =0. Expanding Wj m x = 
we get 



a;i(0,0, Ai) = -a;— h/xoc 



wi(l,0, Ai) = -X 



Ai 
4Ai 



fli 



2 V Ai 



8 Ui 



1-^ ^' + 0{x^ 



which yields at c — > 00 (putting U = c = 1 again) 



wi(0,0,Ai)=/.o 
Analogously, we shall have 

^2(0,0, A2) = -Mo 

a;2(l,0, A2) = -Ho 



,wi(l,0, Ai) = /xo 



0(a;3) 



f gai 



8 V Ai 



(16) 



(17) 



^ _ 1 / g{ai + 02) 
A2 

1 f g jai + a2) 
A2 



^3(0,0, A3) = Ho 



1 f ga2 

A3 



,a;3(l,0, A3) = Ho 



goA' 

8 V A3 y 



(18) 
(19) 



where, of course. A, = ±1 and A^ = 1. 
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Table 1: Gauge coupling constant, mass parameter fiQ and parameters of the 
confining SU(3)-connection for charmonium. 
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Mo «! a2 


bi 


b2 




B2 




(GeV) 


(GeV) 


(GeV) 






4.68010 


0.627818 0.0516520 2.45565 


-0.705320 


1.70660 


2.10247 


4.25862 



3 Relativistic spectrum of charmonium 

Now we can adduce numerical results for constants parametrizing the charmo- 
nium spectrum which are shown in table I. 

One can note that the obtained mass parameter fiQ is consistent with the present- 
day experimental limits |H] where the current mass of c-quark (2/Lto) is accepted 
between 1.1 GeV and 1.4 GeV. As for parameters Ai^2 of solution (5), only the 
wave functions depend on them while the spectrum does not and within the 
present paper we consider Ai = A2 ^ 0. 

With the constants of Table I, the present-day levels of the charmonium spec- 
trum were calculated with the help of (9)-(ll) while their nonrclativistic values 
with the aid of (17)-(19), according to the following combinations (we use the 
notations of levels from ref.[S]): 

rjailS) : si = ui{0, 0, -1) + 1^2(0, 0, -1) + c^3(0, 0, -1) , 

J/i^{lS) : £2 = toi(0, 0, -1) + 6^-2(0, 0, 1) + 603(0, 0, -1) , 
XcoilP) ■■ £3 = wi(0, 0, -1) + co2(0, 0, -1) + ^3(0, 0, 1) , 

Xci(l-P) : £4 = wi(0, 0, 1) + 0^2(0, 0, 1) + 6^-3(0, 0, 1) , 
ijcilP) ■■ £5 - ^1(0, 0, 1) + 6^-2(1, 0, -1) + 0-3(1, 0, -1) , 

Xc2ilP) ■■ £6 = ^1(0, 0, -1) + C02(l, 0, -1) + 6^3(1, 0, -1) , 

V,{2S) : £7 = wi(0, 0, 1) + 6^-2(1, 0, 1) + ^^3(1, 0,-1) , 

^(25) :£8 =wi(0, 0,-1)+ CJ2(1, 0,1) +6^-3(1, 0,-1), 

^-(3770) : £9 = LJi{l, 0, -1) + UJ2{1, 0,-1) + ^^3(0, 0, 1) , 

^-(4040) : £10 = u;i{0, 0, 1) + 6^2(0, 0, -1) + ^^3(1, 0, -1) , 

7/^(4160) : £11 = uji{0, 0, 1) + UJ2{0, 0, 1) -I- ^^3(1, 0, -1) , 

^-(4415) : £12 = uji{0, 0, 1) + ^^2(0, 0, -1) + ^^3(1, 0, 1) . (20) 

Table II contains experimental values of these levels (from ref.[S]) and our the- 
oretical relativistic and nonrclativistic ones, and also the contribution of rela- 
tivistic effects in %, where it makes sense to speak about such a contribution. 
Besides, one can notice that the form of the wave functions (13)-(14) permits to 
consider, for instance, the quantity l//3i to be a characteristic size of quarkonium 
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Table 2: Experimental and theoretical charmonium levels. 



^3 


Experim. 

(GeV) 


Relativ. 
(GeV) 


Nonrelativ. GeV 
(GeV) 


Relativ. contrib. 

(%) 




r 
(fm) 




(fm) 


ro/r 


£l 


2.97980 


2.97980 


2.37202 


20.3965 


0, 


.0603473 


1 


.32755 


21.9984 


£2 


O.U9DOO 


o.uyoo ( 




Z0.4UDU 


U 


.UDUo4 ( 


i 


.OZ ( 00 


zi.yyo4 


£3 


3.41730 


3.41729 


2.37202 


30.5877 


0, 


.0603473 


1 


.32755 


21.9984 


£4 


3.51053 


3.51764 


2.37202 


32.5678 


0, 


.0602972 


1 


.32755 


22.0167 


£5 


3.52614 


3.53085 


1.05011 


70.2590 


0, 


.0602972 


1 


.32755 


22.0167 


£6 


3.55617 


3.54759 


1.05011 


70.3993 


0, 


.0603473 


1 


.32755 


21.9984 


£7 


3.59400 


3.65282 


1.05011 


71.2521 


0, 


.0602972 


1 


.32755 


22.0167 


£8 


3.68600 


3.66955 


1.05011 


71.3831 


0, 


.0603473 


1 


.32755 


21.9984 


£9 


3.76990 


3.75160 


-30.0324 




0, 


.0655399 


2 


.64046 


40.2878 


£10 


4.04000 


4.04906 


33.4683 







.0602972 


1 


.32755 


22.0167 


£11 


4.16000 


4.16614 


33.4683 







.0602972 


1 


.32755 


22.0167 


£12 


4.41500 


4.41872 


33.4683 







.0602972 


1 


.32755 


22.0167 



state. Under the circumstances, if one calculates l//3i in both the relativistic 
(61 7^ 0) and nonrelativistic (5i — 0) cases, then one can obtain those sizes r 
and ro in fm (1 fm = lO^^'^ cm) so the latter are adduced in table II together 
with the quantity r/ro. 

4 Physical interpretation 

The results obtained allow us to draw a number of conclusions. As is seen from 
table II, relativistic values are in good agreement with experimental ones, while 
nonrelativistic ones are not. The contribution of relativistic effects can amount 
to tens per cent and they cannot be considered as small, as was expected by a 
number of theorists [3]. Moreover, the more excited the state of charmonium 
the worse the nonrelativistic approximation. For very excited states, the latter 
is not applicable at all. The physical reason of it is quite clear. Really, we 
have seen in the nonrelativistic limit (see the relations (16)-(19)) that the pa- 
rameters bi,2,Bi 2 (see eq. (5)) of the linear interaction between quarks vanish 
under this limit and the nonrelativistic spectrum is independent of them and 
is practically getting the pure Coulomb one. As a consequence, the picture of 
linear confinement for quarks should be considered as an essentially relativis- 
tic one while the nonrelativistic limit is only a rather crude approximation. In 
fact, as follows from exact solutions of SU(3)-Yang~-Mills equations of (5), the 
linear interaction between quarks is connected with colour magnetic field that 
dies out in the nonrelativistic limit, i.e. for static quarks. Only for the moving 
rapidly enough quarks the above field will appear and generate linear confine- 
ment between them. So the spectrum will depend on both the static Coulomb 
colour electric field and the dynamical colour magnetic field responsible for the 
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linear confinement for quarks which is just confirmed by our considerations. 
In our case, the interaction effect with the colour magnetic field is taken into 
consideration from the very outset just reflects the linear confinement at large 
distances. 

Also, one can notice from table II that r^/r » 1 for all the charmonium states, 
which additionally points out the importance of the relativistic effects connected 
with colour magnetic field for confinement. 

Finally, I would like to say a few words concerning the nonrelativistic potential 
models often used in quarkonium theory. The potentials between quarks here 
are usually modelled by those of harmonic oscillator or of funnel type («. e., of 
the form a/r + /3r with some constants a and (3), see, e. 3., refs. |Hlini- 
It is clear, however, that from the QCD point of view the interaction between 
quarks should be described by the whole SU(3)-connection = A'^Tc, gen- 
uinely relativistic object, the nonrelativistic potential being only some compo- 
nent of A1 surviving in the nonrelativistic limit at c — *■ c». As is easy to show, 
however, the connection of form ~ Br'', where i? is a constant, may be solu- 
tion of the Yang-Mills equations (4) only at 7 = — 1, i. e. in the Coulomb-like 
case. Consequently, the potentials employed in nonrelativistic approaches do 
not obey the Yang-Mills equations. The latter ones are essentially relativistic 
and, as we have seen, the components linear in r of the whole A^ are different 
from At and related with colour magnetic field vanishing in the nonrelativistic 
limit. That is why the nonrelativistic potential approach seems to be inconsis- 
tent though it developed a number of techniques and physical interpretations (in 
particular in charmonium theory (HJ |U] which can be useful under a relativistic 
description as well. Our approach uses only the exact solutions of Yang-Mills 
equations as well as in atomic physics the interaction among particles (e. g., 
the electric Coulomb one) is always the exact solution of the Maxwell equations 
(the particular case of the Yang-Mills equations). 

5 Concluding Remarks 

As we have seen, the application of the Dirac equation to the charmonium 
specrtum leads to a reasonable physical picture. From the obtained results 
there follows that the standard approach of potential models on the base of 
the Schrodinger Equationwith some potential modelling confinement seems to 
be inconsistent. The more consistent approach could be o the basis of the 
Schrodinger Equation in (colour) magnetic field since the linear confinement at 
large distances could be connected with the colour magnetic field rather than 
with the static colour electric one which follows from the exact solution of the 
SU(3)-Yang-Mills equations (see ref. 0] and the solution (5)). Historically, the 
latter way was rejected due to incomprehensible reasons. The analysis of the 
present paper shows, however, that the most consistent approach is probably 
the one based on the Dirac equation in the confining SU(3)-Yang-Mills field 
when the theory is relativistic from the very outset. In its turn, this approach is 
the direct consequence of the relativistic QCD Lagrangian since the mentioned 
Dirac equation is derived just from the latter one. 

The calculations of the present paper can be extended. Indeed we have the 
explicit form (7) for the relativistic wave functions of a quarkonium that may 
be applied to the analysis of the quarkonium radiative decays. Besides, our 
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preliminary calculations show similar results to hold true also for bottomonium. 
At last, there is a possibility of modifying the gluon propagator on the basis of 
exact solutions of the SU(3)-Yang-Mills equations described here and in ref. j3] 
for the mentioned propagator to be able to lead to linear confinement between 
quarks at large ditances. The author hopes to discuss the mentioned questions 
elsewhere. 
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